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ABSTRACT 

We present the first detailed X-ray observations, using Chandra , of NGC 6482 - the near- 
est known 'fossil group'. The group is dominated by an optically luminous giant elliptical 
galaxy and all other known group members are at least two magnitudes fainter The global 
X-ray properties (luminosity, temperature, extent) of NGC 6482 fall within the range of other 
groups, but the detailed properties show interesting differences. We derive the gas temperature 
and total mass profiles for the central 30hjQ kpc 0.1r2oo) using ACIS spatially resolved 
spectroscopy. The unusually high Lx /Lgpt ratio is found to result from a high central gas 
density. The temperature profile shows a continuous decrease outward, dropping to 0.63 of 
its central value at 0.1r2oo- The derived total mass profile is strongly centrally peaked, sug- 
gesting an early formation epoch. These results support a picture in which fossil groups are 
old, giving time for the most massive galaxies to have merged (via the effects of dynamical 
friction) to produce a central giant elliptical galaxy. 

Although the cooling time within 0.1r2oo is less than a Hubble time, no decrease in 
central temperature is detected. The entropy of the system lies toward the low side of the 
distribution seen in poor groups, and it drops all the way into the centre of the system, reaching 
very low values. No isentropic core, such as those predicted in simple preheating models, is 
present. Given the lack of any central temperature drop in the system, it seems unlikely that 
radiative cooling can be invoked to explain this low central entropy. The lack of any signature 
of central cooling is especially striking in a system which appears to be old and relaxed, and to 
have a central cooling time < 10^ years. We find that the centrally peaked temperature profile 
is consistent with a steady-state cooling flow solution with an accretion rate of 2 Mq yr~^, 
given the large P dV work arising from the cuspy mass profile. However, solutions involving 
distributed or non-steady heating cannot be ruled out. 

Key words: galaxies: clusters: general - galaxies: elliptical - galaxies: haloes - intergalactic 
medium - X-ray: galaxies - X-rays: galaxies: clusters 



1 INTRODUCTION 

In a cosmological hierarchy where smaller systems form prior to 
the collapse and virialization of more massive systems, groups form 
before clusters. However, distinguishing old groups among all the 
groups observed today has generally proved impossible. In galaxy 
groups, galaxy mergers can occur efficiently because of their low 
velocity dispersion, similar in some cases to the internal velocity 
dispersions of galaxies. In an old, relatively isolated group with 
little subsequent infall, there may be sufficient time for most mas- 
sive galaxies (but not the low mass galaxies) to lose energy via dy- 
namical friction and merge, producing a system consisting of a gi- 
ant central elliptical galaxy, dwarf galaxies and an extended X-ray 
group halo. This is one of the proposed scenarios for the formation 
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of giant isolated elliptical galaxies and has its origins in the galac- 
tic cannibalism ofHausmai^^Ostntei,(J^97g^. The observation of 
'fossil' systems i|PonmaiTet^ni^9'^ gave th e first direct c lue to 
the above mechanism. Numerical simulations ('Barnes"l989') sug- 
gest that a few billion years is required for compact group mem- 
bers to merge and form a single elliptical galaxy. It is then feasible 
to form such a system in a Hubble time, but if it falls into a cluster 
during this process then the evolutionary history will be obscured. 



Ther e are about a do zen such system s iden- 
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been studied in detail. Based on their sp ace density, an d using a 
precise definition of such a system, Jone s et ^1 i2003h estimate 
that fossil systems represent 8%-20% of all systems with the same 
X-ray luminosity, and are as numerous as poor and rich clusters 
combined. Therefore studying these systems gives us information 
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about the formation of giant elliptical galaxies at redshifts which 
are readily accessible and without the additional complications of 
infall into clusters. 

Jones et al. (200 3) found that fossil systems had an unusually 
high Lx/Lopt ratio compared to normal groups, and suggested 
that the high X-ray luminosity was due to a high density of hot 
gas, reflecting an earl y epoch of formation w hen the Universe was 
denser. Alternatively. IVikhlinin et \\99'f) suggested that such 
systems have an unusually high mass-to-light ratio and thus rep- 
resent partially failed groups with low star formation efficiencies. 
If the distinctive properties of fossil groups arise from an unusu- 
ally early collapse epoch, then they should have strongly peaked 
mass profiles (e.g. Navarro et al. 1995). They are expected to con- 
tain central cool gas, because no merger induced mixing has taken 
place recently. It has also b een suggested that these systems should 
contain lower gas entropy Jjones et al. Il2003h . A given amount of 
energy injected at a high density epoch should result in lower en- 
tropy than if it is injected at a later, low density, epoch. 

The discovery of fossils, and measurements of their X-ray 
properties, has so far been based on ROSAT observations. In this 
paper we report on a detailed Chandra observation of the nearest 
known fossil group. 

This paper is organised as following: Section 2 briefly reviews 
the group properties and describes the data and preparations. The 
results from imaging and spectral X-ray analysis are presented in 
section 3. Section 4 describes the distribution of mass and mass- 
to-light ratio. X-ray scaling relations are discussed in section 5. A 
discussion and concluding remarks are summarised in section 6. 



2 OBSERVATION AND PREPARATION 
2.1 The group 

The NGC 6482 group was selected for Chandra observations via 
a search for nearby fossil groups, based on a sample of isolated 
optically luminous elliptical galaxies. The isolation criteria were 
that all galaxies listed in NED within 0.75 Mpc and 1000 km s~^ of 
the elliptical galaxies were required to be more than 2 magnitudes 
fainter. In addition, the X-ray luminosity, based on the ROSAT All- 
Sky Survey, was required to be > 10"^^ erg s~^, much greater than 
normal ellipticals and meeting the ljones et al. 1 j2003h definition of 
a fossil group. 

The closest galaxy [or group] found which satisfied these 
criteria was NGC 6482 (RA. 17:51:48.83 and Dec. -h23:04:18.9; 
J2000), at a redshift z = 0.0131. This is a giant elliptical galaxy 
and brightest group member. There are 5 confirmed members of 
this group within ~ 60 arcmin of the brightest member. These are 
NGC 6482, UGC 11018, MRK 895, UGC 11024 and CGCG 141- 
014. The second brightest galaxy, UGC 1 1024, has a total (2MASS) 
J magnitude of 11.41, compared with Jtot=9.35 for NGC 6482. 
Thus all the known members are indeed more than 2 magnitudes 
fainter than the brightest galaxy. 

There are more than 20 2MASS near-IR galaxy detections 
within ~ 20 arcmin, the area limited to the locations of the above 
members, with a total J-magnitude ranging from 14.2 to 15.5, but 
there is no independent source to confirm their group membership. 
The velocity dispersion of the group is f»242 km s~^ calculated us- 
ing the heliocentric velocities of individual group members (NED, 
NASA/IPAC Extragalactic Database). 

Analysis of a pointed ROSAT PSPC observation showed that 
NGC 6482 was at the centre of extended X-ray emission. The X- 
ray luminosity was Lx,boi ~ 1-3 x 10^^ erg s~^. An X-ray point 
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Figure 1. ACIS raw image of NGC 6482, centred on S3, in the full energy 
range. X-ray contours from ROSAT observation is overlaid. 

source coinciden t with NGC 6482, probab ly an AGN, hampered a 
detailed analysis. iGoudfrooii et al.Ul994h have noted the presence 
of a radio source and emission-line gas and dust in the nucleus of 
NGC 6482, probably also associated with the active nucleus. The 
derived mass of ionised gas and dust is typical of other elliptical 
galaxies in the study. This study also shows that there is no strong 
dust lane in the brightest group galaxy suggesting the absence of 
any recent merger induced mixing. 

Hq = 70 km s"^ Mpc"^ and fi^ = 0.3 are assumed 
throughout this paper At the redshift of NGC 6482 the luminosity 
and diameter distances are 56.7 Mpc and 55.2 Mpc, respectively 
and 1 arcsec = 0.268 kpc. 

2.2 Chandra observation and data preparation 

NGC 6482 was observed using Chandra ACIS-S on 2002 May 
20. The target was located at the centre of chip 7 (S3), a back- 
illuminated CCD. The total exposure time was over 19 ks. Bad pix- 
els were removed using the supplied bad pixel file. Background 
flares were also removed from the event 2 file leaving us with a 
useful exposure time of just over 18 ks. This study is limited to 
the S3 chip which is entirely illuminated by the X-ray emission 
from the source. A fraction of the X-ray emission is lost in the gaps 
between the chips (Fig.0 and closest part of chip 2 (to the NE) 
contains point sources which appear to be the main source of the 
observed elongation in the ROSAT image. For this analysis, unless 
stated otherwise, the standard CIAO v2.3 routines and tools were 
used. 



2.2.1 Spatial analysis 

Based on our ROSAT analysis, it was suspected that the X-ray 
emission extends beyond the chip under study, S3, and that blank 
sky observations must be used to estimate the background. We use 
those provided in the CAIO calibration data base CALDB. The low- 
est flux regions at the edges of the S3, which are the closest to being 
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Figure 2. X-ray contours from the soft (0.3 - 1 keV) diffuse emission over- 
laid on a 8 X 8 arcmin size DSS optical image. 



'source free', show a 30% higher soft X-ray flux on average com- 
pared to the flux from the same regions at the re-projected blank sky 
observation. At the same time the hard X-ray count rate measured 
for the same regions shows no offset, confirming that the diffuse 
source emission has fully covered the chip. The background esti- 
mated from surrounding chips cannot be directly used since they 
are all front-illuminated. Thus a background image from blank sky 
observations with energies from 0.3-2.0 keV was used for the spa- 
tial analysis. 

Point sources were detected using the CIAO wavelet detec- 
tion algorithm and replaced by their local surrounding mean counts 
before any spatial analysis was carried out. Vignetting and other 
sensitivity variations were corrected for, using an exposure map 
appropriate for the energy band. 



2.2.2 Spectral analysis 

We extracted ACIS spectra in successive circular annuli, in the en- 
ergy range limited to 0.5-2.0 keV, excluding point sources. The 
background was chosen for each annulus separately, to account for 
vignetting, and from the same region on the re-projected blank sky 
observations to account for local variations within the chip. 

T he spectra were fit w ith absorbed hot plasma models and the 
APEC ("S mithetal. 1200 ll) model was found to give the best fit. A 
fixed hydrogen column density of AfjY.ja! = 0.08x10^^ cm^^ was 
included in the model to account for Galactic absorption. An inte- 
grated spectrum was first extracted from a large region, of 3 arcmin 
radius. A mean temperature of 0.66 keV and an unabsorbed flux 
of 0.175 X 10"^^ erg s'^ cm'^ (0.5-2 keV) was found. However, 
the fit was poor with a reduced of about 1.5. This high from 
a single temperature fit is not surprising, since we will see below 
that a strong temperature gradient is present in the X-ray emitting 
plasma. Absorption intrinsic to the source was also required by the 
fit, withavalueof 0.14(±0.02) x 10^^ cm"^. 




Figure 3. Hardness ratio map of NGC6482 on a linear scale. The dark 
region at the centre is the hardest region and the white outer region the 
softest. 



3 RESULTS 

3.1 X-ray morphology and surface brightness 

The diffuse 0.3-2 keV X-ray emission, (see Fig|2}, shows a rather 
relaxed X-ray morphology in the core of the group. For the pro- 
duction of this image, point sources were removed, the background 
subtracted, and the 8x8 binned image adaptively smoothed and 
corrected using the exposure map. The X-ray contours in the cen- 
tral region are circularly symmetric but the outer region shows a 
moderate NE-SW elongation. 



3.1.1 Hardness ratio map 

To study any possible substructure and its significance, a hardness 
ratio map was obtained by dividing the adaptively smoothed dif- 
fuse emission image in a hard band by that in a soft band. The 
soft image contains photons with energies from 400 eV to 900 eV. 
The hard image consisted of photons with energy ranging from 900 
eV to 2 keV. We replaced the point sources by their surrounding 
counts before smoothing the hard band 8x8 binned image using 
the csmooth task in CIAO . The same scale map is then used to 
smooth the soft band image as well as the associated background 
and exposure map images (i.e. all these are smoothed in an identical 
way). Before deriving the hardness ratio map, the corresponding 
backgrounds were subtracted from the soft and hard images. The 
hardness ratio map is limited to the central 3 arcmin radius, Fig.|3| 
due to the limited S/N. The relatively uniform and circularly sym- 
metric distribution of the hardness ratio suggests that the system is 
relatively relaxed and has not experienced a violent or strong per- 
turbation recently, as caused by mergers and starforming activity. 
Small scale features at the centre of the map are found to be at the 
noise level and most likely due to the unresolved point sources or 
inaccurate replacement by local diffuse emission after point source 
removal . We therefore assume that the system is spherically sym- 
metric and is in hydrostatic equilibrium. 
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Figure 4. Surface brightness profile of NGC 6482. The sohd fine represents 
the one-dimensional double /3-model excluding the central point source. 
The gray profile is the best single /3-model fit when the central 5 arcsec 
(1.34 kpc) is masked out. Dotted profile shows the observed PSF profile 
normalised to the surface brightness profile at the centre. 



3.1.2 Gas density profile 

The radial surface brightness profile extracted from the inner 200 
arcsec of the 0.3-2.0 keV image, was fitted using a one-dimensional 
/3-model, 

E(r) = Eo[l + (-)']-''''^'^', (1) 

where ro and 13 are the core radius and index, and Eq the central 
surface brightness. 

Under the assumptions of spherical symmetry and isothermal 
gas, the above surface brightness profile implies a 3-dimensional 
gas number density density distribution given by 

n(r)=no[l + (-)^l-^'^/^ (2) 

Fitting the /3-model given by equation 1 , excluding the central 
point source, we found [3 — 0.53 and core radius, ro — 7.5 arcsec, 
but the reduced was large, 7.5. Excluding 5 arcsec from the cen- 
tral region improved the by factor of 2 but it is still far from a 
good fit to the surface brightness at the central region of the system 
which is crucial in the estimation of other physical quantities such 
as mass and gas density. 

Thus we used a double (3 model where an isothermal /3-model 
is fit separately to the inner and outer regions of the radial surface 
brightness distribution The surface brightness in this case is given 
by 

S^(^) = Eol[l + (^)^]-«''^+^/^ r<i?,„, 

(3) 

EaW = Eo2[l + (^)']-'''^+^/^ r > R,^t 

The surface brightness, and its slope, from the above equations 
must be continuous at the cut radius, Rcut, which itself is chosen to 
minimise x^ ■ In addition, the continuity of the gas density deriva- 
tive must b e met in order to derive a unique mass at Rcut. These 
constraints jPratt & A maud 200Z) were directly applied while fit- 
ting the model in Sherpa. 
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Figure 5. Deprojected APEC model fit to the spectra from 7 successive 
annuli extending to a radius of 150 arcsec. The total reduced is l Ol- 



The results from this fit are roi — 2.95 ± 0.027, ro2 ~ 
26.99 ± 0.22 arcsec and /3i = 0.46 ± 0.002, 02 = 0.59 ± 0.004 
with R cut = 40 arcsec. For c omparison, a ROSAT study of the 
system fe anderson et al. l2003h gave a value /3=0.48, and a poorly 
constrained core radius of ~13 arcsec. Fig.|4|shows the radial fit to 
the surface brightness profile from the double /3-model excluding 
the central point source. The gray profile represents the single /3- 
model fit to the surface brightness profile when the central 5 arcsec 
is excluded. The obvious improvement in the fit is also quantified 
by the reduced x^ of 0.9. 

The central gas density is required to fully describe the dis- 
tribution of the gas distribution. This can be obtained from spectral 
analysis of a reasonably large region around the ce ntre. The method 
is described in the following section and further bv lMaughan et al.1 
12003.) . 

3.2 Spatially resolved spectroscopy 

An absorbed hot plasma APEC model was fit to spectra extracted 
from 7 annuli extending to a radius of 150 arcsec. The annuli were 
chosen to contain 1000 to 1 100 net counts, excluding point sources. 
The spectrum of each annulus was rebinned to contain a minimum 
of 20 counts per bin to increase the S/N. The Ancillary Response 
File (ARF), representing effective area as a function of photon en- 
ergy, was calculated for each annulus since it varies across the chip. 
Due to the extra absorption introduced to the low energy band, be- 
cause of the hydrocarbons built up on the CCDs since the telescope 
launch, an additional correction was needed. This correction was 
applied to each individual ARF. The Redistribution Matrix File 
(RMF), which is used to map the photon energy to pulse height, 
was also obtained for each annulus. Since the X-ray emission cov- 
ers the entire chip, spectra of the background were extracted for 
each annulus from the same region of the blank sky observations. 

The absorbed APEC model was fit to the extracted spectra 
from 7 annuli. APEC gave a better fit to the data compared to other 
hot plasma models such as MEKAL. The innermost region fitted 
with a relatively high reduced x^ , and hence more complex mod- 
els were investigated, as discussed in Section 3.2.2. With the cur- 
rent data we are not able to constrain the abundances even when a 
large region is chosen. Thus we assume a similar abundance for all 
the annuli resulting in a value of Z=0.76±0.28 Z p) in XSPEC's de- 
fault, ANGR, system ( Andres & Grevessell989h . Intrinsic absorp- 
tion was again required in the fits. 
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Figure 6. Electron number density profile. The data points are the avraged 
electron density for each shell measured directly from the X-ray emission 
(section 3.3). The solid curve shows the electron number density obtained 
from the double /3-model fit. Dottted curves correspond to the error in the 
central gas density (see text). 



To correct for the effect of projection, de-projection was per- 
formed in XSPEC using PROJCT model. The reduced from the 
fit to the 7 spectra simultaneously was found to be 1.01 with 248 
degrees of freedom. The PROJCT model fit to the spectra from 
all annuli is shown in Fig.|5| Bolometric luminosities and intrinsic 
absorption of the shells are presented in Table 1. The intrinsic ab- 
sorption is consistent with zero in the outermost annuli, as expected 
if it is associated with the galaxy. 

In order to derive the central gas density, a circular region with 
a radius of r\ = 40 arcsec was chosen and an APEC model was fit 
to the extracted spectrum. The APEC normalisation is defined as: 



^Apec — 



10" 



UeUHdV cm 



(4) 



where Da is the angular distance to the source. With the /9-model 
fitted to the inner region we calculate the integral which is now 
written as the following: 



n^nHdV = 1.17 



nH(r)ATTr {l~cos6)dr. 



(5) 



where tanO = ri/r and assuming that = 1.17n//. Knowing 
the APEC normalisation from the spectral fit, we found the central 
hydrogen nummber density, uhqi = 0.151 ± 0.061 cm""^. The 
value of nHQ2 ~ 0.010 cm~^ (see equations 3) was then obtained 
from the continuity of the density at Rcut ~ 40 arcsec. 

The above calculations are based on the assumption that the 
emissivity of the gas simply scales as density squared. Given the 
steep temperature gradient, this is not obviously justified. However, 
at these temperatures, emissivity is not a strong function of temper- 
ature. In order to estimate the effect of temperature variations on the 
value of electron density, we calculate the electron number density 
directly from the emission of each shell (section 3.3) and compare 
with the model 3D electron number density obtained from the dou- 
ble /3-model fit. Fig|5|shows that the measured electron density is 
in reasonable agreement with that obtained from the /3-model anal- 
ysis. 
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Figure 7. Temperature profile of NGC 6482. Bold symbols show the de- 
projected profile. The light symbols represent the projected temperature 
profile. The curve represents the polytropic fit to the de-projected tempera- 
ture profile with 7 = 1.14. 



3.2.1 Temperature profile 

Fig.0shows the annular and de-projected temperature profiles out 
to 150 arcsec. The temperature profile shows a drop outward reach- 
ing 0.63 of its central value at a mean radius of 29 kpc. It is useful to 
have an analytical temperature profile, which can be used to derive 
the entropy profile. We adopt a polytropic model: 



T{r) = T(0)[1 + (Il)2]-3/2/3(7-i)^ 



(6) 



where rc and /3 are the same parameters as one obtains from a f3- 
model to the surface brightness profile. We fixed (3 at an average 
value of 0.53, from the inner and outer regions, and fitted the poly- 
tropic temperature profile, equation 6, with the core radius and 7 
exponent as free parameters. We found 7 — 1.14 and rc = 12.5 
arcsec. The profile is presented in Fig|7| 

3.2.2 Central annulus and the point source 

The innermost region gave a poor spectral fit. Here we inves- 
tigate two possible causes. The temperature profile shows no 
sign of cooler gas in the central region, contrary to many X- 
ray g r oups of similar mass jHelsdon & Ponmanll200o3 : ISun et al~l 
12003: iMushotzkv et al. Il2003h . Since the point source detection 
routine found a point source at the centre of the galaxy in both low 
and high energy bands, we suspected that it might contain an active 
galactic nucleus (AG N). There is fu rther motivation for this in the 
literature - Goudfrooi i et al. I <I994 found an Ha + [Nil] flux of 
10.8(±1.8) X 10"^* ergcm"^ s"^ in the core of NGC 6482. 

To check whether a central point source could be contami- 
nating our spectral fit in the innermost region, which extends to 
r — 4 arcsec, we excluded the central 2 arcsec, which encompasses 
most of the point spread function of Chandra , and monitored the fit 
quality. We found no significant improvement in the fit compared 
to the original reduced of 1.12. The limited number of counts 
from the central point source does not permit any further analysis 
of a possible AGN component. An upper limit for the X-ray emis- 
sion from the point source is estimated assuming a canonical AGN 
contribution (a power law spectral model with photon index=1.7) 
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Table 1. Results from the PROJCT spectral fit to the seven annuli. 



Annulus 


mean radius 


temperature 


Lx.bol 


intrinsic riu 




(?i7okpc) 


(keV) 


10'*"'^ergs/s 


10^2 j.m~2 


1 


0.54 


0.776 ± 0.037 


0.64 


0.177 ±0.069 


2 


1.88 


0.736 ± 0.032 


0.88 


0.080 ± 0.066 


3 


3.48 


0.696 ± 0.053 


0.62 


0.202 ± 0.073 


4 


5.76 


0.637 ± 0.045 


0.71 


0.164 ±0.066 


5 


9.25 


0.677 ± 0.049 


0.65 


0.119 ±0.075 


6 


14.47 


0.563 ± 0.029 


0.65 


0.003 ± 0.050 


7 


28.94 


0.481 ± 0.035 


1.05 


0.000 ±0.081 



in the spectrum of the innermost region. The spectrum is consistent 
with a maximum AGN luminosity of « 1.0 x 10*" ergs/s within 
90% confidence interval (i.e. /S.^^ =2.71). 

A second possibility is that gas in the high density core is cool- 
ing out. We therefore included a cooling flow component, and fitted 
the central annulus with an absorbed APEC+CFLOW model. This 



fitted with M = 3.21 



h2.50 
-2.06 



M0 yr and a slightly increased in- 



trinsic column of uh ~ 0.19 x lO'^^ cm~^, giving an improved 
reduced of 1 -03 with 30 degrees of freedom. The implications 
of this are discussed later in the paper. 




Figure 8. Entropy profile from 7 successive annuli. Solid and dot-dashed 
lines are the entropy profiles assuming polytropic and isothermal temper- 
ature profiles, respectively. The thick grey line is the mean profile for a 
sample of cool (0.3-1.3 keV) systems derived bv lPonman et al, ]|2003'), and 
dashed line represents the typical profiles seen in cosmological simulations 
including only gravity and shock heating. 



3.3 Entropy 

The entropy of the X-ray emitting gas is defined here as: 

S'(r) = fcT(r)/ne(r)^''^keVcm^ (7) 

where rie is the electron density. The polytropic model described 
above was used to obtain the temperature at a given radius. Sim- 
ilarly the electron density was obtained from the the /3-model gas 
density profile. The resulting parameterised entropy profile, scaled 
by 1/T (where T =0.66 keV is the mean temperature) to facilitate 
comparison with other groups and clusters, is shown in Fig|8|as the 
solid line. 

For comparison with other systems, it is also useful to scale 
such profiles to the virial radius, defined here as the radius, r2oo, 
within which the mean total density of the system is 200 times the 
critical density of the Universe at the present time. In section 4. 1 
below, we will derive the value of r2oo from the present Chandra 
data, by extrapolating fitted models for gas density and temperature 
to derive a total mass profile. However, this involves an extrapola- 
tion of these models by a factor of ~ 10 in radius, beyond the 
extent of the Chandra spectral data which we fit. We therefore ^k- 
fe r (except in section 4. 1) t o use the value r2oo = 361 kpc, derived 
bv'Sande rson et aL ]b003) from ROSAT spectral imaging data ex- 
tending to a radius of 100 kpc, which has been derived excluding a 
large region at the centre, and this has been used to scale the x-axis 
in FigH 

To estimate the average entropy in each shell directly from 
the data, the density was obtained in the following way. The de- 
projected luminosity of each 3D shell was combined with the model 
emissivity e at the shell temperature (and in the same energy band, 
0.3-2.5 keV) and the shell volume V , via: 



L — UennVe, 



(8) 



to derive the electron density. This was then combined with the 
deprojected temperature, to give a spot value for the entropy. The 
resulting entropy values and their associated errors are presented 



in FiglHl and agree well with the analytical curve. The dot-dash 
line shows the entropy profile assuming an isothermal model, and 
clearly overestimates the entropy outside the innermost regions, 
which is reduced by the steep decline in temperature. 

Under the assumption of self-similarity of cluster properties, 
S{r /r2oo) profiles should scale as l/T, so we can compare directly 
with scaled models from simulations of rich clusters. Such simula- 
tions , as well as analytical m odels of spherical accretion onto clus- 
ters jTozzi & NormarJl200ll) . predict entropy profiles with slope 
S oc r^'^, outside a central region in which the profile may flat- 
ten. The dashed line in Fig|H|shows such a profile, normalised us- 
ing simulations by Scott Kay (private communication) including 
only gravitational processes and shock heating. Clearly the entropy 
of the gas in NGC 6482 is substantially higher than would be ex- 
pected from such self-simila r scaling, as seems to be the case in all 
groups jPonman et al. l2003[) . 

The m ean S/T profil e derived for a set of cool (0.3-1.3 keV) 
systems bv lPormian et al. 1 12003) is also shown in Fig[8| and lies 
somewhat above t he profile seen in N GC 6482. However, as ha s 
been pointed out bv lSun et ain<2003h and lMushotzkv et al.ljioosi) . 
a large amount of scatter is seen in the entropy within galaxy 
groups. Comparison with other entropy profiles of poor groups 
based on data from XMM-Newton and Chandra , is instructive. 
Converting where necessary to Ho = 70 km s^^ Mpc^^, the re- 
sult for 5'(0.1r2oo) reported for NGC 1550 (T = 1.37 keV) by 
Sun et al. 1 2003) is 120 keV cm^ whilst M ushotzkv et al. I <2003l) 
find values of 107 keV cm^ and 250 keV cm^ respectively in 
NGC 4325 (T = 0.95 keV) and NGC 2563 (T = 1.36 keV). The 
entropy of NGC 6482 at the same normalised radius is 89 ± 12 keV 
cm^, which is lower than any of the above. 

However NGC 6482 is cooler than these other groups. 
IPoimian et al. I j2003) find that the entropy scaling law in groups 
and clusters is well represented, on average, hy S (X rather 
than the self-similar relation S (x T. Applying a i/yO-^^ scajjjjg 
to the four groups discussed above, we find that NGC 6482 has a 
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Figure 9. Cooling time profile of the core. The dotted line represents 12 
Gyr. 



scaled entropy well below that of NGC2563, but slightly higher 
Xhan NGC 1550 anA NGC4325. 

A striking feature of the profile in Fig|8| compared to the en- 
tropy profiles seen in most simulati ons, or in the averaged obser- 
vational profiles derived by Ponman et al J i2003h . is the absence 
of any significant core to the entropy distribution. The entropy ap- 
pears to drop all the way into the centre, to the limit of the reso- 
lution of our Chandra data, correspondin g to 10~ '^ r2oo. S imilar 
behaviour was reported in NGC 1550 bv lSun et al.1 <2003h . how- 
ever, in NGC 1550 the temperature drops inside 0.05r2oo> so that 
the continuing entropy decline could be due to the effects of radia- 
tive cooling, whilst in NGC 6482 the temperature rises continously 
into the centre. 



3.4 Gas cooling time 

The cooling time of the gas was calculated by dividing the total 
thermal energy of each shell by the associated de-projected lumi- 
nosity. The resulting profile is presented in Fig.|9| The cooling time 
of the entire core, within 0.1r2oo, is less than a Hubble time, and it 
drops to ~ 10* years within the region resolved by Chandra . 




r (Kpc) 

Figure 10. The best fit NFW mass profile is shown by gray (green) sohd 
curve. The concentration parameter is c 60. The dashed-dotted and long 
dashed lines represent the NFW profiles with lower concentrations 40 and 
20, respectively. The tick solid curve shows the total gas mass profile. All 
the profiles extend to 361 kpc. The errors in the mass were estimated using 
MC simulations. The dotted line is the mean mass from the MC simulation. 



iNeumann & Bohringeil il995h have suggested a Monte-Carlo 
method to derive the empirical mass profile and the associated er- 
rors. The gravitational mass consists of two parts. One is the varia- 
tion in gas density, d in n/d In r and the other is that of the temper- 
ature, d\nT / d\nr. For the latter we have generated 1000 phys- 
ical temperature profiles with a temperature window of 0.05 keV 
and step size of 8 arcsec within the observed temperature range. 
A physical temperature is one which guarantees a monotonically 
increasing mass with radius. In this simulation the density profile 
parameters were fixed at the values from the /3-model fit. In order to 
estimate the contribution of the density variation (dlnn/dlnr) to 
the mass error, the surface brightness profile at the location of each 
data point was fitted with d In n/d In r as the free parameter instead 
of 13. The error was estimated at the 68% confidence level using the 
error matrix provided by the fitting program. We checked the re- 
sults by monitoring the variation. In this experiment the values 
of other parameters were optimised. The total error in the mass at 
each data point was then derived by adding the two uncertainties 
quadratically. 



4 GAS AND TOTAL MASS PROFILES 

We use the gas density profile and temperature profile to derive the 
gas and total gravitational mass of the system assuming hydrostatic 
equilibrium and spherical symmetry. The total mass is given by: 



kT{r)r dlnp{r) dlnT{r) 



+ 



dlnr 



(9) 



Gfirup dlnr 

where G and rrip are the gravitational constant and proton mass and 
fi = 0.6. 

It is straightforward to derive an analytical description of the 
mass distribution using our analytical expressions for density and 
temperature. However estimating the error in the mass for the data 
points is not similarly simple, since there are contributions from 
uncertainties in both temperature and density, and the constraints 
of physical validity should ideally be taken into account. 



4.1 NFW mass profile and concentration 

Motivated by the results of numerous c osmological N-body simu- 
lations, we attempt to fit a NFW profile iNavarro et al. Il995h 



4pm(rs) 



(r/r,)(l + r/r,)2 



(10) 



to the total gravitational mass density. We can then obtain the fol- 
lowing integrated mass profile for a spherical mass distribution. 



Mtot{< r) = 1&-Kpm{rs)rl[rs ln(l + r/vs) - 



1 -\- r /rs 



(11) 



where pm{rs) is the density at rs. The mass concentration param- 
eter can then be defined as c = r2oo/''s. By definition, r2oo is the 
radius at which the mean gravitational mass density is 200 times 
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Figure 11. Cumulative mass-to-light ratio profile. The errors are estimated 
from the errors in mass and luminosity for each shell. The solid curve was 
obtained by dividing NFW profile by the luminosity profile of the brightest 
group galaxy. 



the critical density, Pc{z). For our assumed cosmology (Sim, =0.3, 

A=0), the critical density is pc(z) = + 2)^(1 + zfi™). 

We fit the NFW mass profile, equation 11, to the data pre- 
sented in Fie llOl to obtain Ts and derive the the mean density pro- 
file. The best fit profile has =5.1 kpc, and is shown in Fig llOl as 
a thin solid line. r2oo is then calculated by extrapolating the profile 
to 200 pc- This gives r2oo ~ 310 kpc, which is consistent within 
errors with the value of r2oo = SGlljpg kpc derived from ROSAT 
data bv lSanderson et al. and used elsewhere in the present 

paper. The extrapolated total mass at r2oo is M200 ~ 4 x 10^^ M0. 
This best fit profile therefore has a concentration parameter 
310/5.1 = 61, which is remarkably high. For comparison. Fig 



[T0|also shows NFW mass profiles with concentration parameters of 
40 and 20. These are clearly inconsistent with our data. 



4.2 Gas fraction 

The integrated gas mass can be simply derived by integrating the 
gas density (equation 2), and the result is shown in Fig llOl The most 
reliable measure of the gas mass fraction is that within a radius 
of 0.1r2oo, since the temperature profile and hence total gravita- 
tional mass is limited to this radius. We found /gas(0.1r2oo}.= 
0.02. A similar gas fraction is seen in NGC 1550 1 Sun et al~l 
l2003h and NG C 4325, but is sigiiifican tly higher than that found 
in NGC 2563 iMushotzkv et al. I I2OO3I) . Extrapolating implies a 
gas fraction which rises to 0.05 at 0.3r2oo, approximately dou- 
ble the typ ical value derived at this radius for cool systems by 
ISandersonetal.U2003l) . 

If we make a major extrapolation (factor of five for our surface 
brightness data, which primarily determines the gas density, and 
factor 10 for our spectral data, which are required to constrain total 
mass) of our models to r2oo, we derive a total gas mass Mgas[< 
»'20o) ~ 6.7 X IQ^^Mq, and infer a gas fraction fgas{r20o) = 
0.16. The latter is at least as large as that seen in rich clusters, 
and certainly considerably higher than that normally reported for 
galaxy groups. 



4.3 Mass to light ratio profile 

With the high spatial resolution of the Chandra data, it is possible to 
derive the mass-to-light ratio profile in the core of the system. We 
analysed archived R-band CCD imaging data from the 0.9m JKT 
to model the cumulative luminosity profile of the brightest group 
galaxy. The photometric accuracy, based on calibration using five 
nearby stars in the field, was «0.1 mag. 

Fig ll ll shows the cumulative mass-to-light ratio profile within 
0. lr2oo ■ To obtain a continuous profile (the solid curve), we divided 
the fitted NFW mass profile by a numerical model fitted to the R- 
band luminosity of the galaxy. 

The profile shows a rapid rise in mass-to-light ratio with ra- 
dius in the central 5 kpc. Taking into account the luminosity of 
the other known group members, the mass-to-light ratio in the 
R band at r2oo rises to 71 ± 15 Mq/Lq. This is equivalent to 
86±15 Mq/Lq in the B band, assuming colours of a typical el- 
liptical. This value is within the scatter of the mass-to-light ra- 
tios found by Sanderson & Ponman ( 2003) for groups and clusters, 
and is significantl y smaller than the R-band mass-to-light ratios of 
250-450 found bv jVikhlinin et al. 1 1999) for similar systems. The 



IVikhlinin et al. I il999h mass measurements were based partly on 
temperature estimates from an Lx — T relation, and are thus less 
reliable than our mass measurements based on an accurate temper- 
ature profile. 



5 SCALING RELATIONS 

The X-ray scaling relations of groups differ from those of clus- 
ters. Here we investigate whether this fossil group follows the scal- 
ing relation s of g roups, or whether it supports the suggestion of 
I Jones et al. I i200 3). that there are significant differences between 



fossil groups and normal groups. 



5.1 X-ray luminosity 



From the spectral analysis we find the total X-ray luminosity, 0.5- 
8.0 keV, to be 5.21 x 10''^ ergs/sec within ~0.1r2oo- Beyond this 
the total X-ray luminosity of the group has been estimated by ex- 
trapolating the surface brightness profile from the /3-model fit to 



r2oo- We find Lx,i 



1.11 X 10 ergs s within 7-200- 



5.2 The Lx ~ Lopt relation 

The brightest group galaxy has a total R-band luminosity of 4.46 x 
10^" L0, based on our analysis of the archived JKT observation. 
This means that NGC6 482 follows t he sam e trend as the six other 
fossil groups studied in IJones etal.N2003ft . in having an X-ray lu- 
minosity significantly higher than normal groups of the same opti- 
cal luminosity. 

The high X-ray luminosity in NGC 6482 is clearly due to the 
high gas density and is not a result of a high mass-to-light ratio. 
The central density of 0. 15±0.06 cm~^ is higher than that in any 
group or cluster in the compilation of .Sanderson et al. . l2003.) . 



5.3 The Lx -T relation 

Ijones etal.lj2003h found that two of their fossil groups were over- 
luminous compared to the expectations from the cluster Lx — T 
relation, in contrast to the majority of groups which may be un- 
derluminous compared to the cluster relation. With an overall 
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temperature of 0.7 keV the NGC 6482 system falls within the 
scatter of previous measurements of the group Lx-T relation 
jMulchaev & Zabludoff"'l998': Helsdon & Ponman "^OOOb) and is 
not significantly offset. Thus not all fossil groups are offset from the 
Lx-T relation. A possible explanation may be that th e NGC 6482 
system was essentially optically selected, whereas the I Jones et al. I 
y003|) systems were X-ray selected, and their study may therefore 
have preferentially selected high Lx fossil groups. 



6 DISCUSSION AND CONCLUSIONS 

In this section we first give a summary of the results from this de- 
tailed study of the nearest fossil group using high resolution Chan- 
dra data and discuss some of them in detail: 

• The system has extended X-ray emission similar to bright X- 
ray groups but the optical light is dominated by the central giant 
elliptical galaxy. 

• The uniform and circularly symmetric distribution of the hard- 
ness ratio suggests that the system is relatively relaxed with no sign 
of recent violent perturbation. 

• This is the first Chandra resolution observation of a relaxed 
system in which the temperature decreases continually outwards 
despite a central cooling time of only ~ 10* years. 

• In agreement with recent studies of low temperature systems, 
the entropy profile shows no flat core. 

• The Mass-to-light ratio is 71 ± 15 Mq /Lq (R-band) at r2oo. 

• The mass concentration parameter has a remarkably high 
value of c ~ 60, suggesting an early epoch of formation. 

The proximity of NGC 6482, coupled with the high angular res- 
olution of Chandra , has allowed us to resolve the X-ray structure 
of the system in a way never previously achieved with any fossil 
group. We find that the high value of Lx / Lopt, in this fossil group 
at least, results from a very high central gas density (0.15 cm^'^), 
whilst the mass-to-light ratio, measured accurately here, is not ex- 
ceptional. 

The internal absorption may be due to the substanti al gas and 
dust inferred from optical data iGoudfrooii et al. II 1994b . We esti- 
mate a hydrogen column density of ~ 1.7 x 10^" cm~^ based on 
the optical extinction, well below the absorption derived from our 
X-ray spectra. A plausible explanation may be the presence of a dif- 
fuse 'hot dust' component in the central few kpc, as suggested by 
IRAS & ISO data i n many ellipticals ( Goudfrooii & de Jong 1995; 
iFerra ri et al. '2002^. The small size dust grains of such a component 
would be undetected using optical extinction measurements. 

The centrally-concentrated gas density profile, coupled with a 
T(r) profile which rises continuously into the centre of the group, 
implies a remarkably centrally-concentrated total density profile, 
with an NFW concentration parameter c ~ 60. This suggests 
that the system was formed at an early epoch, has been relatively 
undisturbed since, and thus reflects the high density conditions of 
the early Universe. Although N-body simulations involving cold 
dark matter predict that low mass systems such as NGC 6482 (we 
measure a total mass of about 4 x 10^^ Mq) should form before 
rich clusters, and hence have higher concentration parameters, the 
pro perties of NGC 6482 s eem altogether exceptional. For exam- 
ple, 'Wec hsler et al. I <2002 ) identified a total of 14,000 halos in a 
cosmological simulation, and studied their properties. They found 
strong correlations between halo mass, formation epoch and con- 
centration parameter, in the sense that lower mass halos tend to 
form earlier, and have higher c. However, the highest values of c 
for halos from their simulation, in the mass range 3-4 x 10^^ M©, 



was c ~ 25. Such halos formed (according to their definition of 
'formation time', based on the rate of fractional mass growth) at a 
redshift z ~5. 

How are we to understand the extremely high concentration of 
NGC 6482? It is possible that this is pointing to some deficiency 
in cosmological simulations. Some previous authors have reported 
surprisingly high values of c, compared to those expected from 
simulations. For example David et al. (2001) found c ~ 12 in the 
moderately rich cluster A780, which is three times greater than the 
expected value assuming a AC DM cosmology, a nd would require 
an abnormally early redshift of formation of z m 4. |Pratt & Amaud 
(2002) found a mass concentration of c ~ 5.4 in their study of 
A1413 by fitting an NFW profile. Similarly they obtained c ~ 3.75 
for A1983 (.Pratt & Arnaud..2003.) . In the present case, at least, it is 
worth reflecting that the scale radius of 5. 1 kpc lies well within the 
central galaxy, in a region where (c.f. Fig. lilt the stellar mass is 
the dominant component. Hence the scale radius, and the concen- 
tration parameter calculated from it, may be significantly affected 
by dissipative processes taking place in the baryonic component 
of the system. Such processes are known to be poorly modelled in 
cosmological simulations, where they are included at all. For ex- 
ample, it is quite conceivable that inflow of gas into the centre of 
NGC 6482, such as is believed to take place during galaxy mergers, 
could deposit large quantities of cooling gas in the inner regions of 
the galaxy, triggering a nuclear starburst and steepening the central 
mass density profile. Such processes may well have raised the value 
of c subtantially, however the D25 radius of the galaxy is only ~ 16 
kpc, so we can have some confidence that the primordial scale ra- 
dius must have been no larger than this, and therefore that the con- 
centration parameter must have been high (c >20), indicating an 
early formation epoch for the system. 

Thus the NGC 6482 system is an old group, giving sufficient time 
for dynamical friction to act on the most massive galaxies within 
the group, causing them to migrate to the centre, and merge to pro- 
duce the giant elliptical galaxy NGC 6482, leaving no other lumi- 
nous galaxies within 2 mag of the dominant galaxy. 

As we have seen, the entropy of NGC 6482 at r = 0.1r2oo 
appears to be lower than the average for groups, but not out- 
standingly so, considering its low mean temperature. The most 
interesting feature is the way in which the entropy drops all the 
way into the central 1/2 kpc of the galaxy, where our observa- 
tions become limited by the spatial resolution of Chandra . This 
con flicts directly with the expectations of s imple preheating mod- 
els (Evrard & Henrv '1991; 'Kasiser"l99l'; Cavaliere et al. 'l99l 
Balogh et al. 1999; Valageas & Silk 1999; Tozzi & Norman 20o3) 
in which gas is placed on a high adiabat b efore the formation of 
the system, setting a lower limit (argued bv lPonmMi etal.Ul999h 
to be ^ 100 keV cm'^) below which its entropy cannot drop, 
in the absence of significant cooling. Similar results have been 
reported recent ly by a number o f authors (^Ponman et al. 2003; 
IPratt&ArnauS 12003; Sun et al. 1 120031; iMushotzkv e t al. 20o£, 
however, in the case of NGC 6482, the behaviour is all the more 
striking, given that the temperature of the IGM rises continuously 
inward, so that the low central entropy cannot be attributed to the 
effects of cooling. It appears that simple global preheating mod- 
els are no longer tenable, though this does not prohibit models in 
which local heating, either inside or outside forming groups and 
clusters, is responsible for breaking the similarity scaling between 
galaxy systems of different mass. Even NGC 6482 has an entropy 
profile which lies well above that expected from self-similar scaling 
of clusters. 

Despite the fact that the entire core shows a cooling time less 
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than the Hubble time, and in the innermost regions resolved by 
Chandra it drops to ~ 10* years, there is no evidence for any cooler 
gas at the centre of the galaxy. As a result of recent studies with 
XMM-Newton and Chandra tP etere on et al. 2001 : Bohringer et al. I 
[2002; Sakellio iTeral. ll200A IPeterson et al. IEoObI) it is now well- 
established that gas in the cores of most clusters does not appear to 
be cooling at the rates often inferr ed in the past. The r easons for this 
are still a topic of lively debate jpabian et al. 1 *2001). Despite the 
general lack of really cool gas, a decline by up to a factor ~ 2-3 in 
temperature is commonly seen in the cores of clusters and groups. 
The lack of any decline at all in NGC 6482, immediately eliminates 

one suggested mechanism for counteracting central cooling. Ther- 
I I I I 

mal conduction I Voist et al. 2002) cannot act to transport energy 
into the core against the temperature gradient, a nd so is ruled out. 
A hig hly inhomogeneous metallicity distribution iMorris & FabianI 
I2OO3 I) can lead to a lack of lines from cool gas phases, and hence 
limit the apparent temperature drop in the central cooling regions, 
but it cannot suppress th e signature of cooler gas altogether. Heat- 
ing by cluster merging tiKempner et al. I2003t 1 iGomez et al. .l2002h 
can disrupt cooling g as, and inject energy wh ich may suppress 
cooling for some time iKnight & Pormianll997h . however it seems 
most unlikely that this could explain the situation seen in such an 
apparently old, relaxed system as NGC 6482, given the symmetric 
hardness ratio map. 

Three sources of heating remain, which might help to offset the 
observed radiative energy losses 3.5 x lO*'^ erg s~^ from within 
the central 10 kpc). These are PdV work, supemovae and AGN. 
Given the highly peaked nature of the total density profile, gas will 
be subject to considerable P dV work as it flows inward through 
the core of the system. To explore the effects of this, we seek a 
steady state cooling flow solution by solving the energy conserva- 
tion equation (Fabian et al. _1984.) . 



1 .2 



pv — {H{r) + 0(r)) = n,{r)nH{r)e{T), 



(12) 



where H(r-) 



5 fcT(r) 
2 /xm/j 



is the specific enthalpy and (f) is the gravi- 



tational potential, coupled with mass conservation 

M — A-Kpvr^ = constant 



(13) 



for the observed gas density (equation 2) and the gravitational po- 
tential inferred from our data, to calculate the temperature profile. 

We fix the gas properties at the outer edge of our analysis (where 
the cooling time is approximately equal to the Hubble time) to 
match our observations, and seek for a value of M which produces 
a r(r) profile matching that which we observe (Fig. ll2t . using the 
integral equation 



kT{r) 



Id 



l.nnH{x)e(T)x^dx\^C,(\A) 



where ry — 28.94 kpc is the mean radius of the outermost shell 
and the constant of integration C is chosen to normalise the tem- 
perature profiles to the value observed at r = r-j. A good match to 
our profile is obtained (Fig. ll2> with M ~ 2 M0 yr^ ^ . Lower mass 
accretion rates, M « 1.5, result in a nearly isothermal core, whilst 
higher rates give too steep a profile. Under this steady state cooling 
flow model, the gas is heated by the gravitational work done on it 
as it flows inward, before cooling and dropping out at the centre 
of the galaxy. The spectral results shown discussed above in Sec- 
tion 3.2.2, do not disallow such a cooling component in the centre 
of the system - in fact, as we have seen, it actually improves the 
spectral fits. 

Steady-state cooling flow models have fallen out of favour in 
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Figure 12. Model temperature using the observed mass and gas density pro- 
files for different mass accretion rates. The data point are the de-projected 
temperatures in FiglTl 



recent years, as discussed above, and it is, of course, possible that 
some processes other than P dV work are responsible for maintain- 
ing the steep temperature profile. Any steady-state solution other 
than that which we have derived above would have to involve dis- 
tributed heating, since the gas must be prevented from cooling 
throughout the region we have analysed. 

An obvious source of heating is supernovae (SN), which could 
be of type I and/or II. To estimate the contribution of SNIa, we use 
a total SN I rate of (0.18 ± 0.05)/i?5 per century per 10^° Lbp, and 
an en ergy production of lO^'^ erg per supernovae iCappellaro et al. I 
Il999ft . This corresponds to Lsn/Lb = 2.5 x 10^" erg s"^ Lg^ 
after correcting for the value of Hq used in this study, giving an 
SNIa contribution of ^ 1.0 x 10*^ erg s"\ Hence SNIa could 
contribute a significant fraction of the ~ 3 x 10'*^ erg s"'^ radiated 
by the region under study. 

Using the Ha luminosity taoudfrooii etal.lll994h and associ- 



ated star formation rate, SFR(M(7)/vr~^)=7.9 x 10"*^ l^Ha (erg 
s~^) jCardiel et aLll2003h . we estimate the SNII contrib ution as- 
suming a SN rate 9.26 x 10~^ times the star formation rate iKawatal 
2001), and an energy of 10^ erg per SN. This corresponds to a heat 
injection of 0.14 x lO''^ erg s~^, which is much smaller than the 
power required to replace the radiative losses. 

The final source of heating which might play a role is injec- 
tion of energy from an AGN. AGN heating is achievable either 
by Compton heating, mechanical heating or viscous dissipation 
of pressure wav e s gen erated by AGN (Ruszkowski et al. 2003). 
iBinnev & Taboil il995') have argued that the last heating mech- 
anism is most efficient. These authors, and several others, have 
shown that even a small fraction of AGN mechanical energy trans- 
ferred to the ISM is sufficient to balance the radiative loss of energy 
via X-ray radiation. The bubbles and jets released from AGN are re- 
sponsible for distributed heating, but it is important to note that in 
reality such processes take place in a non-steady fashion. 

A key feature of our observations is the negative temperature gra- 
dient with radius. Heated cooling flow models are able to reproduce 
such a negative temperature gradient (e.g. Brighenti & Mathew^ 
J2OO2?), but a detailed comparison is required before more defini- 
tive conclusions can be drawn. Thus a non-steady state solution 
remains a possibility. However for a steady state solution, we have 
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shown that PdV work alone could generate a temperature profile 
similar to that observed in NGC 6482. 
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